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ABSTRACT 

^' 
O , 

O , We report late-time observations of the radio afterglow of GRB 030329. The 

light curves show a clear achromatic flattening at 50 days after the explosion. 
We interpret this flattening as resulting from the blast wave becoming trans- 
relativistic. Modeling of this transition enables us to make estimates of the energy 
i content of the burst, regardless of the initial jet structure or the distribution of 

initial Lorentz factors of the ejecta. We find, in accordance with other events, 
^ i that GRB 030329 is well-described by an explosion with total energy of a few 

O i times 10^^ erg expanding into a circumburst medium with a density of order 

O 

00 

o 

Q ! Subject headings: gamma-ray bursts: specific (GRB 030329) - radio continuum: 

^ ! general - cosmology: observations 

o 

c/3 . 1. Introduction 

To fully understand the progenitors of gamma-ray bursts (GRB) and the workings of 
^ I the central engines that power them, we must first understand their energetics. The wide 
I dispersion in isotropic GRB energies, once corrected for their jet-like geometry, results in a 
narrowly-clustered radiated energy ~ 1.3 x 10^^ erg (Bloom, Frail & Kulkarni 2003; Frail 
et al. 2001). Likewise, the kinetic energy Ek of the outflow, as derived from the afterglow 
emission (Wijers & Galama 1999; Kumar 2000; Freedman & Waxman 2001), is clustered 
around a similar value (Panaitescu & Kumar 2002; Berger, Kulkarni & Frail 2003). The 
close agreement between the values of E^ and E^ have lead to the prevailing view that GRB 
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explosions draw from a standard energy reservoir with a total budget of order a few times 
10^^ erg. 

GRB 030329 appeared to challenge this view. A distinct break in the optical and X-ray 
light curves at i ~ 0.5 d (Price et al. 2003; Tiengo et al. 2003) if interpreted as a jet, gives 
an opening angle 9j — 5° and a radiated energy (corrected for the narrow opening angle) 
— 5x 10^^ erg - nearly 30 times below the "standard" value. GRB 030329 is also an order 
of magnitude under-energetic as measured by the geometrically-corrected X-ray luminosity 
(Berger, Kulkarni & Frail 2003). Other low energy outliers in the distribution (Bloom, 
Frail & Kulkarni 2003) and X-ray luminosity distribution do exist. These energy outliers are 
important for understanding the diversity of cosmic explosions. They may represent genuine 
low energy events, a long-lived central engine, the viewing of a jet slightly off axis, or events 
for which the efficiency of converting the kinetic energy in the shock to radiation is low. 

For GRB 030329 two different mechanisms have been proposed, both of which increase 
increase Ek to about 10^^ erg. To explain the continued rise of the radio emission over the 
following week, Berger et al. (2003) introduced a wide-angle jet {9j = 17°) with a lower 
Lorentz factor than the narrow jet but carrying ~ 10 x more energy. This led Berger et 
al. (2003) to conclude that GRBs share a common energy scale only when the total energy, 
including that in mildly relativistic ejecta, is considered. Another variant of the standard 
energy hypothesis is to include a long-lived central engine. This was first proposed to account 
for an extended bright phase of afterglow emission of GRB 021004 (Fox et al. 2003). Granot, 
Nakar & Piran (2003) and Piran, Nakar & Granot (2004) have argued that such model could 
also explain the bumps in the optical light curves (Lipkin et al. 2004) and the rising radio 
light curves (quahtatively at least). 

The above discussion shows that there are grounds to question the simplest hypothesis 
for GRB explosions. Further progress on the nature of the central engines will most certainly 
involve measuring and modeling the temporal and angular behavior of the energy release. 
Nonetheless, the total energy yield of these explosions is perhaps the most important bulk 
parameter of cosmic explosions. Fortunately, as has been noted in the past (Frail, Waxman 
& Kulkarni 2000), late-time radio observations offer us precisely this opportunity. Here, 
we report on observations of the radio afterglow of GRB 030329 on a timescale of 60 days 
to 366 days. Using these observations in conjunction with the earlier radio data (Berger 
et al. 2003), we are able to infer a geometry- independent estimate of the total energy of 
GRB 030329. 
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2. Observations and Results 

All observations were made with the Very Large Array (VLA)^ and reduced following 
standard practice (e.g., Frail, Waxman & Kulkarni 2000). A summary of the flux density 
measurements are in Table 1. In Figure 1 we plot the radio light curves supplemented with 
our earlier data from Berger et al. (2003) and Sheth et al. (2003). 

A significant flattening of the late-time radio light curves at t ^ 50 d can be seen, 
especially when compared to the steep decline measured at early times at higher frequencies 
(22.5-250 GHz) where the temporal slope a ~ -2 (with F oc T) at i > 10 days. This 
flattening is achromatic, except at 1.43 GHz where the light curve was still rising, reaching a 
maximum near 300 d. A joint temporal and spectral power-law fit of the form F(i/, t) oc f^v^^ 
for i ^ 75 d at 4.86, 8.46, 15 and 22.5 GHz gives a = -1.20 ± 0.06 and (3 = -0.58 ± 0.08. A 
similar flattening of the X-ray emission is seen by XMM-Newton in which the decay index 
changed from ox = —1.86 ± 0.06 to ax = —1.40 ± 0.15 between 1.24 and 37 days and 
37 and 258 days after the burst (Tiengo et al. 2003; Tiengo et al. 2004). The late-time 
X-ray spectral index f3x = 1-17 ± 0.04. The difference in the radio/X-ray spectral slopes, 
/3 — Px = 0.59 ± 0.09 is of the right magnitude {Af3 = 0.5) if these bands lie on either side 
of the synchrotron cooling frequency, i.e. ur < Uc < I'x- 

The simplest explanation for an achromatic flattening of afterglow light curves at late 
times is that the blast wave has become trans- relativistic. Such a transition occurs when 
the rest mass energy swept up by the expanding shock becomes comparable to the initial 
kinetic energy of the ejecta. Given the jet parameters of GRB 030329, this is predicted to 
occur roughly at a time tjvj? ^ ^j/^j — 60 d (Frail, Waxman & Kulkarni 2000). 

In the non-relativistic regime the temporal index in a constant density medium (ISM) 
is ajvR = (21 — 15p)/10 for i/ < i>c or a^R = (4 — 3p)/2 for v > here p is the energy index 
of the electrons accelerated into a shock, and it is related to the spectral slope p — —2P for 
1/ > i/c, and p — 1 — 2(3 ior u < Uc- In contrast, in the relativistic regime, the flux is expected 
to steeply decline, t'^ independent of frequency. Similar expressions exist for a a wind-blown 
medium (WIND) whose density varies as the inverse square of the radius: a^R — {5 — 7p)/6 
for u < Uc or unr — {8 — 7p)/6 for u > Uc (Livio & Waxman 2000). 

Applying the above expressions for v < Vcto the radio measurements of a and /3 gives 
p = 2.16 ± 0.16, aNR = -1.14 ± 0.24 (ISM) and aNR = -1-69 ± 0.19 (WIND). Thus the 
observed flattening of the radio light curves {a — —1.20 ±0.06) is consistent with a constant 
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density medium, but a wind-blown medium predicts a decay that steeper than expected. 
A wind-blown medium is rejected at the 2.5a level. Applying the same expressions to the 
X-ray measurements where u > Uc, we derive p = 2.34 ± 0.08, aNR = —1-51 ± 0.18 (ISM) 
and a^TK = —1.40 ± 0.09 (WIND). By themselves, the X-ray data is not able to distinguish 
between models while the radio favors the ISM model. 

Given this good agreement (both in the timescalc and the magnitude of the flattening) 
with the expectations of a jet expanding in a constant density medium, we are motivated to 
consider more detailed models in the following sections. 

3. Modeling the Broadband Afterglow 

Initially the blast wave is expanding relativistically and thus we can use the extensive 
machinery developed for afterglow modeling (e.g., Wijers & Galama 1999; Chevalier & Li 
2000; Panaitcscu & Kumar 2002). The details of our specific method are given in Yost et al. 
(2003). Briefly, the blast wave is assumed to accelerate a power-law distribution of electrons 
N(7e) oc 7"'' with energy index p above some minimum energy 7^. The emitted spectrum 
(synchrotron and possibly inverse Compton) at any time t can be described by the peak flux 
and several break frequencies whose evolution depends on the shock dynamics. 

Prior to the jet break time tj the evolution is identical to a spherically symmetric shock 
(Sari, Piran & Narayan 1998), while for tj < t < t^R the evolution is described by a 
jet expanding in a constant density medium (Sari, Piran & Halpern 1999). After tNR the 
evolution of the break frequencies is governed by Sedov- Taylor dynamics (e.g.. Frail, Waxman 
& Kulkarni 2000) and discussed in the next section. To allow a smooth transition of the 
light curves, the normalization of the emission spectrum is determined by the equations for 
the jet model ai t — tNR- 

A total of six physical parameters were fitted for; the isotropic equivalent energy £'52 
of the jet, its half opening angle 9j, the circumburst density n, the electron energy index p, 
and the fraction of the shock energy in electrons ee and in magnetic fields e^. We used only 
the radio and X-ray data since the optical data are dominated by the supernova 2003dh at 
late times. 

We show the results of our fit in Figure 1 and summarize the model parameters in Table 
2. Also in Table 2 is the best-fit jet model from Berger et al. (2003). The calculation of 
6j from Berger et al. (2003) was originally derived using the isotropic gamma-ray energy 
and some assumed radiative cfliciency. To be internally consistent we have re-computed 6j 
in Table 2 directly from the E^2 and n using the formalism of Sari, Piran & Halpern (1999). 
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4. Modeling the Afterglow at Late Times 

As noted above, after tNR the blast wave becomes non-relativistic. On about the same 
timescale the side- ways expansion of the jet becomes important and eventually the blast 
wave becomes spherical. Modeling the afterglow emission at late stages thus offer the distinct 
advantage of the inferred parameters being independent of the relative orientation of the axis 
of the jet and the hne of sight to the observer. 

We follow the methodology of Frail, Waxman & Kulkarni (2000) (Appendix A) and 
Berger et al. (2004). Essentially the primary observable is the synchrotron self-absorbed 
radio spectrum with the peak flux density (fu), the optical depth (tj/) and the characteristic 
synchrotron frequency (i^m)- These parameters, in turn, are determined by the values of 
four physical parameters at t — ^nr^ the magnetic field (Snr), the shock radius (rNR), the 
minimum Lorentz factor of the electrons (7m,NR) and the number density of the radiating 
electrons (rie.NR)- 

The evolution of the self-absorbed spectrum is governed by the dynamics of the blast 
wave. In particular, at late times the blast wave evolution is well described by the Sedov-von 
Neumann- Taylor (SNT) solution (Zel'dovich & Raizer 2002): r = TNR(t/iNR)^^^; speed (in 
units of speed of light, c) p = /9NR(tANR)"^^^ B = 5nr(^Anr)"^''^ and 7^ = 7m,NR(^AjvR)"^/' 
We drop the NR subscript hereafter for clarity. 

There are four unknown parameters but usually only three observables (but see §4.3). 
Thus we must introduce additional constraints to fully specify a solution. We consider three 
possible constraints and these are discussed in turn. 



4.1. Equipartition 

The standard approach in radio astronomy is to appeal to equipartition for an additional 

constraint. This means equating the fraction of energy in relativistic electrons. Eg = [{p — 
l)/(p — 2)] X Ue^rti'iTT'eC^V to the energy in magnetic fields, Eb = B^V/Sti or eg = e^. Here, 
V is the volume within the shell of emitting electrons, V = Anr^/r] cm^, where r/r] is the 
thickness of the shell and we use the standard (thin shell) assumption that 7] ~ 10. The 
density n = He x 3/r] where 3/?] is the shock compression factor. 

Under the assumption of equipartition we find the following best-fit solution for rest 
frame physical parameters: r^p ~ 3.5 x 10^^ ^m%^^^^ cm, implying B^p ~ 0.18 z^^V^^ ^> 
jm,ep ~ 42 z/O^^g^ ne,ep ^ 5.9 1/"°/^ cm'^ and €e,ep i^^^g = €B,ep i^^^g ~ 0.08 for values of 
p ^ 2.2 and ^nr ~ 50 days. 
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As a consistency check, we note that applying Equation 13 of (Frail, Waxman & Kulkarni 
2000) gives f3ep ~ 1.25 i^^^g^^^^ and thus we confirm that it is appropriate to apply the SNT 
solution at ^nr- 

Because i/m drops below our lowest frequency just before t-^R, the model parameters 
are given in terms of i'm,9 which is normalized to 1 GHz at ^nr for p — 2.2, in accordance 
with our relativistic modeling (§3). Using the model constraints above we derive £'min = 
i?e,ep + Es^ep ~ 1.3 X 10^° ^mV^^ ^^S- The total energy in the ejecta (including shocked 
protons), derived from the Sedov- Taylor energy equation is £^sT,ep ~ 8.5 x 10^° ^^g^^ erg. 

4.2. Consistency with the Sedov- Taylor Solution 

Instead of adopting the equipartition assumption in §4.1, we can use the more general 

constraint that the energy in the electrons and magnetic field does not exceed the thermal 
energy (Frail, Waxman & Kulkarni 2000), + Eb < i^sx/S. This constraint is a logical 
consequence of our assumption that the blast wave can be described by a Sedov- Taylor 
solution at late times. 

For radii above (below) r^p, the summed energy rises steeply as r^^ (r~^) when Eb > E^, 
{Eb < Ef.). With this method we can probe the parameter space at(;a?//rom equipartition and 
derive a range of allowed values for Est- Figure 2 shows the dependence of -Est and [E^+Eb) 
on radius. As can be seen from the Figure, the permitted values of the Sedov- Taylor energy 
are E'st ~ 0.64 — 2.2 x 10^^ erg, corresponding to radius values of r 2.2 — 4.0 x 10^^ cm. 
This range of radii imphes the following ranges for the physical parameters: B 0.027 — 0.30 
G, rie ^ 2.3 - 160 cm-^, (3 ^ 0.78 - 1.4, ee ^ x 10"^ - 0.49, ^ 1.7 x 10"^ - 0.45. 

4.3. The Fourth Observable - CooUng Frequency 

Late-time X-ray observations (Tiengo et al. 2004) as in §3 allow us to constrain the 
cooling frequency z/^, and hence determine the radius of the blastwave. In particular. X-ray 
observations allow us to model the evolution of the cooling frequency, Uc oc B'H"'^ (Frail, 
Waxman & Kulkarni 2000). For the same value of v^^g used above: r «i 3.6 x 10^^ cm, 
B ^ 0.21 G, -fm « 39, He ^ 4.3 cm-^ and €e ^ 0.064, ^ 0.13 (Table 2). The total energy 
in the ejecta Est ~ 7.8^fg x 10^° erg. The uncertainty in the energy is dominated by the 
fiux density error of the X-ray observations at t ~ 60 days. 
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5. Angulcir Size of the Afterglow 

Comparison of our derived size with the observed size offers the opportunity to directly 
verify models. Normally such a comparison is inaccurate since the expected angular sizes are 
small and until recently, they were inferred from interstellar scintillation observations (e.g., 
Frail et al. 1997). However, Taylor et al. (2004) have directly determined the size of the 
fireball for GRB 030329, measuring an angular diameter 9 — 172 ± 43 //as at t — 83.3 d. 

Using Sedov- Taylor expansion (i.e., r oc t^/^) to scale to tNR=50 d, and assuming a 
Lambda cosmology {Hq = 71 km/s/Mpc, Q,m = 0.27, il\ = 0.73), the angular-diameter 
distance of GRB 030329 at ^ = 0.1685 is dA = 589 Mpc and its radius r = (6.2 ± 1.5) x 10^^ 
cm. Thus it appears that the VLBA measurement is at odds with the radius estimates 
deduced above, albeit at only 2a significance (see Figure 2). To better understand this 
difference (and to extract more robust estimates of r and EsT/n), it would be useful to 
jointly fit the evolution of the apparent size of GRB 030329 (Oren, Nakar & Piran 2004) 
together with the evolution of the afterglow emission. 

6. Discussion and Conclusions 

Continued monitoring of the afterglow of GRB 030329 has revealed an achromatic flat- 
tening (§2) of the light curves whose magnitude and timescale are consistent with a jet-like 
outflow which is undergoing a dynamical transition to trans-relativistic expansion. We have 
modeled the entire dataset (1-366 d) and derive a total kinetic energy of 9.0 x 10^'^ erg (§3). 
This estimate agrees with a previous estimate of 6.7 x 10^° erg, based on observations from 
1-64 d (Berger et al. 2003). Thus our analysis is consistent with earlier conclusions that 
GRB 030329 was not an under-energetic burst (Berger et al. 2003; Granot, Nakar & Piran 
2003). 

We also attempted to derive an independent estimate of the energy using only those ob- 
servations beyond t > tNR- The Sedov-Taylor analysis (§4) traces all ejecta, regardless of the 
initial jet structure or the distribution of initial Lorentz factors of the ejecta. This accurate 
calorimetric method has now been applied to GRB 030329, GRB 970508 and GRB 980703 
(Frail, Waxman &: Kulkarni 2000; Berger et al. 2004). In this particular case, the value 
of this approach is limited by the uncertainty in the exact value of the synchrotron peak 
frequency Um- Using the X-ray data as an added constraint (§4.3), and adopting Um—l GHz 
at t]\fR we derive an energy of 7.8 x 10^° erg. 

Strictly speaking the kinetic energies given above are upper limits since they are derived 
at a flxed time {tNR)- The flreball undergoes radiative losses especially prior to the time 
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when the evolution is in the fast coohng regime (Sari, Piran & Narayan 1998). We estimate 
from our full model in §3 that this fast-to-slow cooling transition occurs at t = 0.03 d. At 
this time the total isotropic energy E^2 = 2.0 x 10^^ erg and the geometrically-corrected 
energy = 2.2 x 10^^ erg. The energy (isotropic or otherwise) drops to 42% of this value 
at tj and 37% at Inr- 

With further observations of GRB 030329 we have an opportunity to look for departures 
in the standard blastwave model. At present, the physical parameters of the blastwave 
derived from only the first two months (Table 2 and Bcrger et al. 2003) provide a remarkably 
good prediction of its subsequent behavior for up to a year later. Deviations may be expected 
to occur at late times because of an injection of additional energy by slower moving ejecta, 
an evolution in the microphysical parameters of the shock (p, eg, e^) with time, or a change 
in the density structure of circumburst medium. 

The current estimates of the energy and other physical parameters of GRB 030329 are 
accurate at best to a factor of two. While a more quantitative analysis (e.g., Yost et al. 2003) 
is needed, we can already rule out order of magnitude changes in E^, Ce a-nd e^. Identifying 
changes in the density structure is more difficult because between 1 day and 1 year the blast 
wave samples only a small range of radii in its rest frame (~0.1 to 0.4 pc). However, at its 
present rate of decay the afterglow of GRB 030329 will be visible at centimeter wavelengths 
for the next decade. Thus the prospects for tracing the continued evolution of GRB 030329 
are excellent. 

DAF wishes to thank Pawan Kumar for his hospitality at the University of Texas, 
Austin during which time this paper was begun. We acknowledge NSF and NASA grants 
for support. 
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Table 1. Observational Summary 



Date 


At 


-^1.43 


-^4.86 


-^8.46 


-Pis.o 


-^22.5 


(UT) 


(days) 


(/^Jy) 


(/^Jy) 


(/^Jy) 


(/^Jy) 


(/^Jy) 


(1) 


(2) 


(3) 


(1) 


(5) 


(6) 


(7) 



2003 Jun. 10.00 
2003 Jun. 15.00 
2003 Jun. 25.01 
2003 Jun. 30.94 
2003 Jul. 09.02 
2003 Jul. 27.88 
2003 Aug. 06.06 
2003 Aug. 11.96 
2003 Aug. 25.89 
2003 Sep. 18.79 
2003 Sep. 22.74 
2003 Sep. 27.75 
2003 Oct. 06.73 
2003 Oct. 30.72 
2003 Dec. 02.58 

2003 Dec. 07.58 

2004 Jan. 31.53 
2004 Feb. 05.36 
2004 Feb. 26.27 
2004 Mar. 06.33 
2004 Mar. 29.41 



72.52 


1009±116 


77.52 


1111±84 


87.53 




93.46 


1360±40 


101.54 




120.40 




129.58 


1276±69 

-L ^ f \J —1—\J tJ 


135.48 




149.41 


1677±66 


173.31 


3470±65 


177.26 




182.27 


2240±66 


191.25 


2450±58 


215.24 


1709±52 


248.10 


1276±40 


253.10 


1541±41 


308.05 


1240±54 


312.88 


1800±53 


333.79 


1360±96 


342.85 


1075±124 


365.93 


904±102 



4808±67 3454±69 

4613±66 3572±54 

• ■ • 3016±31 
3290±57 2130±55 

• ■ • 2299±65 

■ • ■ 1190±35 
1955±62 

■ • ■ 1525±56 
1677±55 1277±51 
1590±52 1122±45 

• ■ • 1002±91 
1600±46 1050±38 
1407±44 931±41 
1154±43 874±44 

999±32 637±33 

949±32 696±31 
686±38 

1029±26 675±24 

865±33 547±27 

886±45 447±33 

767±34 528±29 



2328±207 
2680±200 

2100±165 1650±65 

1729±182 1290±65 

997±171 794±119 

436±200 374±166 



*Table columns are (1) Epoch of observation; (2) time in days from burst, (3) flux density 
at 1.43 GHz, (4) flux density at 4.86 GHz, (5) flux density at 8.46 GHz, (6) flux density at 
15 GHz, and (7) flux density at 22.5 GHz. 



Note. — 



To maximize sensitivity, the full VLA bandwidth (100 MHz) was recorded 
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in each of two hands of circular polarization. The amplitude scale was calibrated 
with observations of J1331+303 (3C286), while phase calibration was achieved with 
frequent observations of J1051+2119 at all frequencies except for 1.43 GHz where 
J1016+2037 was used instead. 
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Table 2. Best Fit Parameters For Different Models 





Range 




tj 


tNR 




9, 


n 








Model^ 


(d) 




(d) 


(d) 


(erg) 


(rad) 


(cm~^) 


P 




(%) 


Rel. 


1-64 


5133/164 


10 


N/A 


0.67 


0.45 


3.0 


2.2 


0.19 


4.2 


Full 


1-366 


3245/241 


14 


48 


0.90 


0.47 


2.2 


2.2 


0.17 


7.4 


Non-Rel 


50-366 


5313/75 


N/A 


50 


0.78 


N/A 


1.3 


2.2 


0.06 


13 



^From top to bottom, the parameters are given for the relativistic model by Berger 
et al. (2003), the full evolution model from §3, and the non-relativistic model from 
§4.3. 

''The large for these fits is affected by intensity variations induced by density 
inhomogeneities in the ionized interstellar medium of our galaxy (Goodman 1997). 
Diffractive scintillation dominates at early times, while the fluctuations in the late- 
time 1.43 GHz hght curve are likely due to large scale refractive effects. 

'^The blastwave kinetic energy in units of 10^^ ergs. The first two models have been 
corrected for coUimation by multiplying the isotropic equivalent energy by (1 — cos^j). 
The energies for the first two models are calculated at tj, while the last model is the 
Sedov-Taylor energy calculated at tNR- 
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Time After Burst (days) 

Fig. 1. — Radio light curves of the afterglow of GRB 030329. From bottom to top the 
measurements are were taken at 1.43 GHz (green square), 4.86 GHz (blue triangle), 8.46 
GHz (red circle), 15.3 GHz (black diamond), 22.5 GHz (green diamond), 43 GHz (blue 
circle), 100 GHz (red triangle), 250 GHz (black square). The data have been scaled relative 
to the 8.46 GHz measurements by u^, where f3 is the spectral slope between 4.86 and 22.5 
GHz for t ^ 75 d (§2). The solid lines are the predicted emission from a jet expanding into 
a constant density medium (see §3). The flattening of the light curves beyond t~50 days 
occurs when the jet expansion becomes sub-relativistic. Although this model captures the 
gross evolution of the observed light curves (the rise to maximum, the initial steep decay 
at high frequencies and the flattening of the optically thin light curves at late times), the 
unprecedented high signal-to-noise of this data does, however, show some deficiencies of 
our model. In addition to the expected fluctuations at low frequencies due to interstellar 
scintillation, there are real variations at higher frequencies such as the "bump" between 50 
and 60 day seen at 15 GHz. Owing to our simplistic treatment, the model also under predicts 
(or over predicts) the flux density at different transition times when the evolution changes 
abruptly. 
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Fig. 2. — Dependence of energy on radius at t = ^nr ~ 50 days as described in Section 4. 
The sum of the energy in electrons and magnetic fields, Eg + Eb, is minimized near the 
equipartition radius. The Sedov-Taylor energy scales with the unknown radius as Est oc r~^. 
The permitted range of kinetic energy values are defined where {Ef. + Eb) < -Est/S and are 
marked by the bounded sohd fine. Here we are adopting i^^^nr =1 GHz. The radius of the 
shock, as measured from VLBI observations of Taylor et al. (2004) at t ^ 83.3 days (and 
scaled down to t ~ 50 days as r oc t'''^) is shown by the vertical shaded band; it is a factor of 
~ 2 times larger (l.So" significance) than that predicted by our NR modeling. Including the 
additional constraint of the X-ray observations (Tiengo et al. 2004), we arc able to constrain 
the radius to 3.6 ± 0.5 x 10^'' cm (vertical shaded band) - nearly equal to the equipartition 
value. This implies a slightly lower energy of Est ~ ^.Sl^ig x lO^'^ erg. 



